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We study CP violation in the two-body decay of a scalar tau into a neutralino and a tau, which 
should be probed at the LHC and ILC. From the normal tau polarization, a CP asymmetry is defined 
which is sensitive to the CP phases of the trilinear scalar coupling parameter A T , the gaugino mass 
parameter Mi, and the higgsino mass parameter fj, in the stau-neutralino sector of the Minimal 
Supersymmetric Standard Model. Asymmetries of more than 70% are obtained in scenarios with 
strong stau mixing. As a result, detectable CP asymmetries in stau decays at the LHC are found, 
motivating further detailed experimental studies for probing the SUSY CP phases. 



I. INTRODUCTION 

The surplus of matter over anti-matter within the 
universe can only be explained with a thorough under- 
standing of CP violation. The CP phase in the quark 
mixing matrix of the Standard Model, which has been 
confirmed by B-mcson experiments [l[, is not sufficient 
to understand the baryon asymmetry of the Universe ■ 
However, the Minimal Supersymmetric Standard Model 
(MSSM) provides new physical phases that are 
manifestly CP-sensitive. After absorbing non-physical 
phases, we chose the complex parameters to be the 
higgsino mass parameter /i, the U(l), and SU(3) gaugino 
mass parameters Mi, and M3, and the trilinear scalar 
coupling parameters Af of the third generation sfcrmions 
(/ = b,t,r). The corresponding phases violate CP and 
are generally constrained by experimental bounds on 
electric dipole moments (EDMs) Q. However, these 
restrictions are strongly model dependent 0-0], such 
that additional measurements outside the low energy 
EDM sector are required. 

Many CP observables have been proposed and studied 
in order to measure CP violation. Total cross sections Q , 
masses 0, and branching ratios are CP-even quan- 
tities. For a direct evidence of CP violation, however, 
CP-odd (T-odd) observables are required. Examples are 
rate asymmetries of either branching ratios [ll|, cross 
sections [lH, or angular distributions [l3j]. Since these 
rate asymmetries require the presence of absorptive 
phases, they are typically small, of the order of < 10%, 
if they are not resonantly enhanced [1J] . Larger CP-odd 
observables which already appear at tree-level arc 
desirable. These are T-odd triple products of momenta 
and/or spins, from which CP-odd asymmetries can be 
constructed. Such triple product asymmetries are highly 
CP-sensitive, and have been intensively studied both at 
lepton and hadron colliders [TH, [lj| . 




the CP phases of the trilinear coupling parameters Af 
arc rather unconstrained by the EDMs 
phases of Ab and A t have been studied in stop 
and sbottom J2{| [26| decays, respectively. Since these 
are decays of a scalar particle, the spin-spin correlations 
have to be taken into account. The triple product 
asymmetries can then be up to 40%, for sizable squark 
mixing. Similarly for probing the CP-violating phase of 
A T in the stau vertex, f~x -T, it is essential to include 
the tau spin. Only then is there a sensitivity to the 
phase of A T [2?], [28( . If the spin of the tau is summed 
over, this crucial information is lost. Triple product 
asymmetries including the tau polarization have been 
studied in neutralino decays \1 ^ tt [28| , and also in 



chargino decays % 



\. It was shown that the 



normal tau polarization itself is CP-sensitive, and that 
the asymmetries are large and of the order of 60% to 70%. 

We are thus motivated to study CP violation, including 
the tau polarization, in the two-body decay of a stau 

f m ^T + xl m = l,2, i = 2,3,4, (1) 

followed by the subsequent chain of two-body decays 

h+L; (2a) 



4^X1 + ^2; n = L,R, e = e,fi. 



(2b) 



See Fig.[T]for a schematic picture of the entire stau decay. 
This process is kincmatically open for a mass hierarchy 



(3) 



where the staus are heavier than the smuons and selec- 
trons. We thus work in MSSM scenarios with heavier 
stau soft SUSY breaking parameters 



M E T > M E C 



M E„ 



(4) 
(5) 



Third generation sfcrmions have a rich phenomenology 
at high energy colliders like the LHC [l7| or ILC [TU due 
to a sizable mixing of left and right states. In addition, 



We show that the normal tau polarization, with respect 
to the plane spanned by the r and l\ momentum, is 
a triple product asymmetry which is sensitive to the 
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FIG. 1: Schematic picture of stau decay. 



phases of A Tl Mi, and /j, in the stau-neutralino sector. 
For nearly degenerate stau masses, Mg ps , a 

strong stau mixing is obtained which results in tau 
polarization asymmetries of more than 70%. This should 
be measurable at colliders 1 . Since the stau is a scalar 
particle, its particular production does not contribute 
to CP-sensitive spin-spin correlations, and can thus be 
considered separately. This allows a collider-independent 
study, where we only discuss the boost dependence of 
the CP asymmetries. 



The paper is organized as follows. In Section |TT] we 
review stau mixing and the stau-neutralino Lagrangian 
with complex couplings. We calculate the amplitude 
squared for the entire stau decay in the spin-density ma- 
trix formalism [30j . We construct the CP asymmetry 
from the normal tau polarization, and discuss its MSSM 
parameter dependence, as well its boost dependence for 
colliders like the ILC and LHC. In Section IIII1 we nu- 
merically study the phase and parameter dependence of 
the asymmetry, and the stau and neutralino branching 
ratios. We comment on the impact of the fi decay in 
scenarios with nearly degenerate stau masses. We sum- 
marize and conclude in Section HVl The Appendices con- 
tain the definitions of momenta and spin vectors, the 
analytical expressions for the stau decay amplitudes in 
the spin-density matrix formalism, and formulae for the 
stau decay widths. 



II. FORMALISM 



A. Stau mixing 



In the complex MSSM, the stau mixing matrix in the 
(f L ,77i)-basis is HEI 



Mf 



-i<j) 



*m T \k?\ 



(6) 



CP violation is parameterized by the physical phase 



4> f = arg[A f ], 
Kf = A T — ji* cot / 



(7) 
(8) 



with the complex trilinear scalar coupling parameter A T , 
the complex higgisino mass parameter fj,, and tan ft = 
V1/V2, the ratio of the vacuum expectation values of the 
two neutral Higgs fields. The left and right stau masses 
are 



- M'^ + (-- + sin 2 9 w )m% cos(2/3) + m 2 , (9) 



1 
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Mp - sin 2 e w m 2 z cos(2/3) + to 2 , 



(10) 



with the real soft SUSY breaking parameters M~ - , 
the electroweak mixing angle 9 W , and the masses of the 
Z boson mz, and of the tau lepton, m T . 

In the mass basis, the stop eigenstates are [!, |3l| 



TV 



TL 



111) 



smf 



with the diagonalization matrix 
TV = 

and the stau mixing angle 
cos Of 



sm Vf 
e~ 1 ^ cos( 



-m T |Af| 



> 2 |A?| 



ml 1 



sm ( 



m 2 fL - m\ 



The stau mass eigenvalues are 



(12) 

(13) 
(14) 



Note that we do not include the tau decay in our calculations. 
However, some of the decay products of the tau have to be re- 
constructed in order to measure the tau spin. The main goal of 
our work is to motivate such an experimental study, to address 
the feasibility of measuring the CP phases at the LHC or ILC. 



n A,2 



(m 2 f L +m 2 f R ) t 



■Am 2 f\kf\ 



(15) 
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B. Lagrangian and complex couplings 



The relevant Lagrangian terms for the stau decay 
-> tx° are d [Hf 

J? Tfi o =gf K„ P R + b f ml P L ) x" f m + h.c., (16) 

with Pl.r = (1 =F 7s)/2, and the weak coupling constant 
g = e/ sin^, e > 0. The couplings arc dchned as [3l[ 



with the ncutralino hclicities A,, A^. The amplitude 
squared decomposes into the remnants of the propaga- 
tors 



A(j) 



Sj - mj + irrijTj 



(27) 



The stau diagonalization matrix 72.* is given in Eq. ([2]). 
and 




BJ = 




(18) 



In the photino, zino, higgsino basis (7, Z, H®, H°), we 
have 

ft = ^ 



sin- ^ N i2 + sin ^iVji 



, (19) 



cos9 w \2 

f% = V2sme w (t a ne w N* 2 ~N* 1 ), (20) 
h L Tl = (ft*)* = -Y T (N? 3 cos/? + AT* sin/3), (21) 



mw cos /3 



(22) 



with mvi/ the mass of the W boson, and N the complex, 
unitary 4x4 matrix that diagonalizes the neutralino mass 
matrix Q 



N* ■ ■ — diag(m^o, . . . , m^o). 



(23) 



The interaction Lagrangian relevant for the neutralino 
decay x° -> ?%. L eT ' for e = e > A* is 

^0 = g£fl i P R x ih+9if£PLX% + h.c. > (24) 
with the couplings /«.' given in Eqs. (|19|) and (|20|) . 



with mass TOj , and width P, of particle j = x® or £, and 
the unnormalized spin density matrices for stau decay 
Pd(j), and neutralino decay PD 1 {x.i)- The decay matrix 
of the spinless slepton is a factor since the polarizations 
of the final lepton and LSP are not accessible. The cor- 
responding amplitude is denoted by D^t). Defining a 
set of spin basis vectors s" for the tau, see Eqs. (|A10[) 
in Appendix [XI and s-p for the neutralino [32j, the spin 
density matrices can be expanded in terms of the Pauli 
matrices a 



Mf)%$ =D6 x ^r6 XiK + E D (<r a ) KK 5 x , K + 

E D 5 X ^ (a») XtX , + E£ fc K) At A; (* 6 )a,A{ , (28) 



PDAXi) 1 = Ui d « + L Di (cr ) ■ 



(29) 



with an implicit sum over a,b — 1,2,3, respectively. 
The real expansion coefficients D, Di, E^, Eq, E d 
and E D h contain the physical information of the process. 
D denotes the unpolarized part of the amplitude for 
stau decay f m — > Xi T 1 Di denotes the unpolarized 
part for neutralino decay Xi ~> ^R^i , respectively. Eq 
gives the tau polarization, E D , and E D describe the 
contributions from the neutralino polarization, and E D b 
is the spin-spin correlation term, which contains the 
CP-sensitive parts. We give the expansion coefficients 
explicitly in Appendix O 



C. Tau spin density matrix 

The unnormalized, 2x2, hermitian, r spin density 
matrix for stau decay, Eqs. ([1]) and ©, reads 



,A T A' _ 



(\M\ 2 ) XtK dJfips, 



(25) 



with the amplitude A4, and the Lorcntz invariant phase 
space element dJzfips, for details see Appendix [Bj The r 
hclicities are denoted by A T and \' T . In the spin density 
matrix formalism [3(| , the amplitude squared is given by 

(|^| 2 ) AtA; = |A(x°)| 2 |A(£)| 2 x 

Y,P°( f t£ p Dl (x°i) KXi D 2 (i), (26) 

A; A' 



Inserting the density matrices, Eqs. (|28|) and ([29]). into 
Eq. (|2"6"|) , we get for the amplitude squared 



{\M\') 



2\a t a; 



(DDi 



,S Dl )<5 



A T A' 



(ESDi+E^E^JK)^^ D 3j (30) 



with an implicit sum over a,b = 1,2,3. The amplitude 
squared (|.M| 2 ) AtAt is now decomposed into an unpolar- 
ized part (first summand), and into the part for the tau 
polarization (second summand), in Eq. (|30j). By using 
the completeness relations for the neutralino spin vectors, 
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Eq. (|A12|) . the products in Eq. (f30|) can be written 2 



D. Normal tau polarization and CP asymmetry 



^Di 



y^afr K^b 
^Di 



mtqipr - Ply) - (Px° ■ Pt)(P£i ■ Px°) , (31) 



{ ^\{\aU 2 + \bU 2 )\fii\ 2 rn T x 



{S°r ■ Py» ) (P X ° ■ Pi! ) - m|o (S" • Pi, ) 

( +/Maf m (&L)*}l/!| 2 "^ X 

(Pt • Px? ) ( S r • Ki ) - (Pr • P*i ) (S? • P*° ) 

+ )5 4 |/#| 2 m x o x 
Jm{a^(6^)*}bf, W 1S Pr, <]• (32) 



The spin-spin correlation term E^ Ej^ , Eq. (|3"2"j) . explic- 
itly depends on the imaginary part 3m{aJ ni (5J ni )*} of the 
stau-tau-neutralino couplings, Eq. (|16[) . Thus this term 
is manifestly CP-sensitive, i.e., it depends on the phases 
4>At, 0i, 4>{i °f the stau-tau-neutralino sector. The imagi- 
nary part is multiplied by the totally anti-symmetric (ep- 
silon) product, 

£ a = \p f , Pll , p T , s a T ] = e^ptpl P p s^, (33) 

with the convention £0123 = 1- Since each of the spatial 
components of the four-momenta p, or the spin vectors 
s", changes sign under a time transformation, t — > —t, 
the epsilon product £ a is T-odd. In the stau rest frame, 
p~ = (mf,0), the epsilon product reduces to the T-odd 
triple product T a 

[Pf, Pei, Pr, s"] = m f (p tl x p T ) • s" = m f T a . (34) 

The task in the next section is to define an observable 
that projects out from the amplitude squared the part 
proportional to £ a (or T a ), in order to probe the CP- 
sensitive coupling combination 3m{aJ ni (b^ ni )*} . 



2 The formulas arc given for the decay of a negatively charged stau 
■-cO Ml™^ h,, ,-o £+l~. The signs in parentheses 



decay t£' -> 



_X?, followed by \" 
in Eqs. J 3 11) and (132 P hold for the charge conjugated stau decay 
t-l&R- ^ n or der to obtain the terms for the 
r ^Xi > however, followed by the neutralino decay 
into a positively charged slepton, ~ ► ^i~^k> one nas *° reverse 
the signs of Eqs. (1 3 1 1 > and J 32 1) . This is due to the sign change 
of Sj^, see Eqs. iC6l . In Appendix [Cl we also give the terms 

for the neutralino decay into a left slepton, jjj? — ► i\t^- Note 
that the term proportional to m r in Eq. I I32H is negligible at high 
particle energies E S> m T . 



The r polarization is given by the expectation value of 
the Pauli matrices cr— (a±, ct 2 , er 3 ) (33j 



(35) 



Tr{p<x} 
Tr{p} ' 



with the t spin density matrix as given in Eq. (|25|) . 
In our convention for the polarization vector "P 
= {Vi, 7*2,1*3), the components 7*\ and P3 are the 
transverse and longitudinal polarizations in the plane 
spanned by pi t and p T , respectively, and 7*2 is the 
polarization normal to that plane. See our definition of 
the tau spin basis vectors s° in Appendix [Al 



The normal r polarization is equivalcntly defined as 
^ _ N(-t)-N(i) 



N(t) + N(l)' 



(36) 



with the number of events ./V with the r spin up (")") or 
down (J,), with respect to the quantization axis pe ± x p r , 
see Eq. (|A10[) . The normal r polarization can thus also 
be regarded as an asymmetry 



V 2 



(j(T > 0) - a(T < 0) 



a(T > 0) + ct(T < 0) ' 
of the triple product 

r = (p ei x Pr ) -£ T , 



(37) 



(38) 



where £ T is the direction of the r spin vector for each 
event. The triple product T is included in the spin-spin 
correlation term E£ 6 E^, Eq. (J32J) , cf. Eq. (JSJ), and 
the asymmetry thus probes the term which contains the 
CP-sensitive coupling combination 3m{al ni {b 7 mi )*} . 

Since under naive time reversal, t — > —t, the 
triple product T changes sign, the tau polarization 7*2, 
Eq. d37J), is T-odd. Due to CPT invariance [H, V 2 would 
thus be CP-odd at tree level. In general, V2 also has con- 
tributions from absorptive phases, e.g. from intermediate 
s-state resonances or final-state interactions, which do 
not signal CP violation. Although such absorptive con- 
tributions are a higher order effect, and thus expected to 
be small, they can be eliminated in the true CP asym- 
metry [ll] 



A™ =~-{V2-V2), 



1 



(39) 



where 7*2 is the normal tau polarization for the charged 
conjugated process f* n — > t + x®. For our analysis at tree 
level, where no absorptive phases are present, we find 
7*2 = —7*2, see the sign change in Eqs. ([31) and (j32j) . 
and thus A^ P = 7*2- Wc study A^ p in the following, 
which is, however, equivalent to 7*2 at tree level. 
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Inserting now the explicit form of the density matrix p, 
Eq. P5|). into Eq. ([31)]), together with Eq. O, we obtain 
the CP asymmetry 



A 



CP 



□ 2 ' b E Dl d^ips 



JDD X dJfips 



(40) 



where we have used the narrow width approxima- 
tion for the propagators in the phase space element 
dJtfips, see Eq. (|D9|) . Note that in the denominator 
of , Eq. (|40|) . the spin correlation terms vanish, 
/Ej) E^dJzfips = 0, see Eq. ((311), when integrated 
over phase space. In the numerator only the spin-spin 



correlation term E^ 6 



E^ for a = 2 contributes, which 



contains the T-odd epsilon product 5 a , see Eq. ([33 




0.2 0.4 0.6 0.8 1 

0f 



E. Parameter dependence of the CP asymmetry 



To qualitatively understand the dependence of the 
asymmetry A^ F , Eq. (|4"0"]) , on the MSSM parameters, 
we study in some detail its dependence on the T m -T-x^ 
couplings, a T mi and b r mi , see Eq. (|D3|) . From the explicit 
form of the decay terms E^ Eq Eq. (I3T1) . and D, Di, 
Eqs. (|C1|) . (|C5p . respectively, we find that the asymme- 
try 



A 



CP 



| d^fips 



(41) 



with (px° -Pt) = (m£ - m| )/2, and (p^o -p^) 



(m| ? 



m~)/2, is proportional to the decay coupling factor 



FIG. 2: Boost distributions of the r polarization asymmetry 
At* ', Eq. ([39]). normalized by A?* (P? = 0), for three different 
sets of stau masses, m T1 2 ~ 200 GeV (solid, red), 400 GeV 
(dashed, green), and 1000 TeV (dotted, blue), see text, for 
stau decay ?i — > rx2, followed by x° - > ^l^fi, and — >• x?^2 
(I = e or p), see Fig. [1] The SUSY parameters are given in 
Table U 



for m — 1, and the sign in parentheses holds for to = 2. 
Thus we expect a maximal rj and thus maximal asymme- 
tries for maximal stau mixing 3 , Of w ±7r/4, and a max- 
imal CP phase in the stau mixing matrix, <f>f rj ±7r/2. 
Note that, in particular, the dependence of <f>f on <f>A T is 
strong for \ A T \ > \ tan/3. We will study numerically the 
phase and parameter dependence on At F and r\ further 
in Section Hm 



Jffl{ g mi( b m«)*} 

|(KJ 2 + fel 2 ) : 



(42) 



with r) mi G [—1,1]. We thus expect maximal asym- 
metries for equal moduli of left and right couplings, 
\ a mi\ ~ \Kni\i which have a phase difference of about n/2, 
where the coupling factor can be maximal r\ m i = ±1, sec 
Eq. (HH). 

To study the dependence of rj on the CP phase <f>f of 
the stau sector, and the stau mixing angle Of , we expand 
the imaginary part of the product of f m -T-x- couplings 

M<i(«4n = 3m{i^ ml i 2 / T L ,^ + i^ m2 | 2 /;f^ 

+ n ml n* m2 [(h^ 2 -f^]}, (43) 

in terms of the stau mixing matrix 1Z, the gauge couplings 

icular, for t 
0, we have 



f T ( R and the higgs couplings h^l ■ In particular, for a 
CP-conserving ncutralino sector, <pi = <j)^ 



3m{aL(&L)*} = Asin0 f sin(20 f U [{h*? - f?J^ 



(-) 



(44) 



F. Boost dependence 

The triple product asymmetry A^ p , Eq. (j40|) . is not 
Lorentz invariant but depends on the boost of the decay- 
ing stau, 



iPfl 



(45) 



In Fig. 2, we show the boost dependence of the asym- 
metry A^ p (/3f), normalized by A^ p (f3f = 0). The 
SUSY parameters are given in Table [I] and we have 
chosen three sets of different f soft-breaking param- 
eters {M£ T ,M ir } = {195,200} GeV (solid, red); 
{395,400} GeV (dashed, green); and {998,1000} GeV 
(dotted, blue). The corresponding stau masses are 



3 Note that a maximal mixing is naturally achieved for nearly 
degenerate staus. However then the asymmetries for fi and T2 
decay typically have similar magnitude but opposite sign, and 
thus might cancel. See the discussion at the end of the numerics 
in Section ITTTDl 
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TABLE I: Benchmark scenario. The mass parameters M2, 
\fi\, A T , Mg, M L M & , and M L are given in GeV. 



Pi 




<pA T 


M 2 


ImI 


At 


tan ft 








tt/2 


250 


250 


2000 


3 










m m 




M- L 


495 




500 




150 




200 



{m fl ,m f2 } = {194,209}; {395,404}; {998,1002} GcV, 
respectively. The corresponding asymmetries in the stau 
rest frame are A^ p (/3?) = -66%; -72%, -71%. Note 
that we have chosen nearly degenerate stau masses which 
lead to an enhanced stau mixing and thus to maximal 
asymmetries; see also the discussion in Section [Oil 



For the stau masses {r 



J = {194,209} GeV, the 



staus can be produced at the ILC with y/s = 500 GeV, 
and have a fixed boost of /3f = 0.63. The corresponding 
asymmetry is then reduced to A^ = —53% if the stau 
rest frame cannot be reconstructed. Typical ILC cross 
section for these masses are of the order of some 20 fb [391 ■ 



If the staus are produced at the LHC, they will have 
a distinct boost distribution depending on their mass, 
which typically peaks at high values (3 T sa 0.9 for stau 
masses of the order of a few 100 GeV up to a 1 TeV, 
see e.g. Refs [2l], [26[. Then the normal tau polarization 
in the laboratory frame is obtained by folding the boost 
dependent polarization A^ P with the normalized stau 
boost distribution (2lj . 



A 



CP 



1 



lab 



dap 



op J d(3f 



A^{fc)dfc, 



(46) 



with the production cross section ap = a(pp —> f + f~). 
The typical reduction of the normal tau polarization 
-^T P iab i s °f the order of two thirds of the asymmetry 
compared to that in the stau rest frame ^ p (0). How- 
ever, it has been recently shown (for similar asymmetries 
in stop decays at the LHC), that the rest frame can be 
partly reconstructed event by event using on-shell mass 
conditions, see Refs. (22|. The LHC cross section for 
stau pair production, cr(pp — > T-j~), also sensitively de- 
pends on the stau masses, e.g., for our benchmark sce- 
nario in Table HI we find cross sections up to 10 fb at 
s/s = U TeV [35|. 



III. NUMERICAL RESULTS 

We quantitatively study the tau polarization asym- 
metry, and the branching ratios for the two-body decay 
chain 



for £ — e,fi. The asymmetry probes the MSSM phases 
(pi, (pfj, and <Pa t i °f the neutralino and stau sector. We 
center our numerical discussion around a general MSSM 
benchmark scenario, see Table |U We choose heavier soft 
breaking parameters in the stau sector than in the e, jl 
sector, to enable the mass hierarchy 



m X.t tR Xi 



(48) 



Further we choose almost degenerate staus which en- 
hances their mixing, leading to maximal asymmetries. 
We choose a large value of the trilinear scalar coupling 
parameter, \A T \ > tan /3 4 , to enhance the impact of 
<PAt in the stau sector. Finally, to reduce the number of 
MSSM parameters, we use the (GUT inspired) relation 
\Mi\ = 5/3 M 2 tan 2 9 W [|[ for the gaugino mass parame- 
ters. The resulting masses of the staus, neutralinos and 
charginos are summarized in Table iHl 



A. Phase dependence 

For the benchmark scenario given in Table U we study 
the phase dependence of the asymmetry A^ in the stau 
rest frame. In Fig. 3(a)| we show the dependence on 
the CP phases in the neutralino sector, <pi and (p^. In 
Fig. 3(b) we show the dependence on the phases in the 
stau secotor <pA T and (p^. The asymmetry strongly de- 
pends on cpA T ~ <pf, which we expect for \A T \ 3> \n\ tan/3 
as in our benchmark scenario, see Table [U In particular 
for (p u = in Fig. [3(b) the asymmetry follows the ap- 

H 



= in Fig. 
proximation formu 
values at d> T 



r> A T 



a Eq 
w ±tt/2 



and attains its maximal 



B. tan/3 dependence and stau mixing 



In Fig. 4(a) we show the \A T \ and tan/3 dependence 
of the asymmetry A^ P in the stau rest frame. We 
can observe that the asymmetry obtains its maximum, 



TABLE II: Mass spectrum for the scenario in Table H] 



m [GeV] 



m [GeV] 



GL, fi>L 
V e , Vn 

V T 

fl 

T2 



155 
204 
192 
497 
495 
504 



X? 

xi 
xf 



112 
190 

254 
327 
181 
325 



X'l'+^V (47) 



4 The value of \A T \ is restricted bz the vacuum stability contidtion 
as \A 2 T < 3(m? + m? x + M% + ti 2 \ [H. 
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37T/2 




37T/2 




FIG. 3: Phase dependence of (a) the r polarization asymmetry At P , Eq. (|39[) . in percent, in the 



plane (for 4>a t = 0), 



and (b) in the 4>a t -4>i± plane (for <j)\ =0), in the stau rest frame. We consider the decay fx -> Txl followed by X l -> 



and ir 



Xi£ 2 where I = e or /i, c/. Fig. [T] The other MSSM parameters are defined in Table [I] 



A^ F ~ —77%, where al so th e coupling factor is maxi- 
mal, r] sa 0.95, see Fig. |4(b)| As discussed in Subsec- 
tion III El the imaginary part of the product of the stau 
couplings 3m{cil ni (tf ni )*} is maximal for a maximal CP 
phase 4>f = 7r/2 in the stau sector, which we show in 
Fig. 4(c) Note that the location of the maximum of A^ F 
is not at maximal stau mixing, sin(0f) = 1/V2 w 0.7, 
since r\ oc sin(2#f )/(|a r | 2 + \b T \ 2 ) starts to decrease for 
increasing A T and tan/3. 

To study the stau mixing, we show the Mg. -Mi de- 



pendence of the asymmetry A^ m Fig- 5(a) In the 



entire Mg, -M^ plane, the CP phase in the stau sector 
is almost maximal, <f>f — 0.617T. However, the asymme- 
try obtains its maxima in the small corridor Mg ps , 
where the stau mixing is maximal, 6f — tt/4. 



C. | /j. | — M2 dependence and branching ratios 

We show the |/x|-M2 dependence of the asymmetry 
A^ F in Fig. [5(b)] The maxima of A^ F are obtained 
where the coupling factor r\ is also maximal, see Eq. (|4"2"j). 

In Fig. [6(a)| we show the corresponding stau branch- 
ing ratio, BR(fi — > TX2), which can be as large as 40%. 
Other competing channels can reach BR(fi 



TX°l) 



65%, and BR(fi -> ^ T xf(2)) ~ 20(10)%. The stau decay 
into the chargino xf is always open since typically the 
second lightest neutralino and the lightest chargino are 



ing ratio BR(^2 —> Hr), summed over I = e, fj,, is shown 
in Fig. |6(b)| which reaches up to 100%. The other im- 
portant competing decay channels are BR(x§ — ^ vivi), 
and BR(x° — > ^l), which open around fj, pa 250 GeV 
and (i sa 300 GeV, respectively, for M 2 = 250 GeV. Note 
that in our benchmark scenario, see Table HI we have 
BR(£ R -> X \t) = I- 



D. Impact of T2 decay 



almost degenerate, m^o 



X-l 



The neutralino branch- 



As we discussed in Section fill B[ we find large asym- 
metries for nearly degenerate staus, where we naturally 
obtain a maximal stau mixing. However, then typically 
the asymmetries for f\ and T2 decay are similar in mag- 
nitude, but opposite in sign. For example in our bench- 
mark scenario we find A^ F = —71% for fi decay, but 
A^ F = +32% for the decay of f 2 . If the production and 
decay process of f\ cannot be experimentally disentan- 
gled from that of f 2 properly, the two asymmetries might 
cancel. We show their sum in Fig. 7(a) in the Mg -M^ 
plane. In Fig. |7(b)| we show the corresponding stau mass 
splitting. 

Note that also the stau branching ratios are similar 
in size; for example in our benchmark scenario we have 
BR(fi -> rxt) = 18%, and BR(f 2 -> txT) = 30%. For 
the Mg -Mi plane shown in Fig. [5j the decay branching 
ratio BR(fi — > rx°) is at least 10%, and that of f 2 is 
larger by roughly a factor of 2 to 4. 
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FIG. 4: |A T |-tan/3 dependence of (a) the r polarization asymmetry _4!j? p , Eq. (|39|) . in percent, in the stau rest frame (for 
the decay fi —¥ ryQ, followed by \2 — > ^J"^, and — ¥ Xi^z tor I = e or fi, cf. Fig. [TJ, (b) the coupling factor r], Eq. ((42}, 
(c) the phase 4>f in the stau sector, Eq. (0, and (d) sin(2#f), with 8f the stau mixing angle, Eqs. I|13|l. (|14|) . The plots are for 
= 7r/4, the other MSSM parameters are given in Table U 
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FIG. 5: Dependence of the r polarization asymmetry , Eq. (|39|l , in percent, in the stau rest frame (for the decay ft — ► rx§i 
followed by ^2 — >• ^x^R) an< ^ ~ * Xi@2 f° r ^ — e or Mi see Fig- Ui on ( a ) the s °ft breaking parameters in the stau sector Mjj, , 
Mi , Eqs. ©, Eqs. (|10l) . In (b) the dependence of A^ F on the gaugino and higgsino parameters \fj,\, M%. Below the contour 
mg R = m-jo the two-body decay X2 — > Ur is kinematically forbidden, above the contour mg R = m^o the lightest neutralino is 
no longer the LSP since mg R < m%o. Below the contour "i^± = 100 GeV the lightest chargino is lighter than 100 GeV. The 
MSSM parameters are given in Tabled] 
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FIG. 6: Contour lines in the \fi\— M% plane of (a) the stau branching ratio BR(ri — y rx®) in percent, and (b) the neutralino 
branching ratio BR(x§ — ► ££r), in percent, summed over both lepton flavors I = e, fj, and charges, for the MSSM parameters 
as given in Table [T] Below the contours nig R = m-o in Figs. 6(a) 6(b) the two-body decay X2 
forbidden, above the contours mg R = m 



the lightest neutralino is no longer the LSP since mg F 



m~± = 100 GeV the lightest chargino is lighter than 100 GeV. 
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FIG. 7: Contour lines of (a) the sum of the r polarization 
asymmetries A^ ¥ , Eq. I)39p . in percent, for the decays fi — > 
r^2 an d ^2 — » t~X2 5 each in the stau rest frame and followed 

by x°2 -»■ <h -» fOT < 

(b) the stau mass splitting mf 2 — ni Tl 
are shown in the plane of the soft breaking parameters of the 



e or n, see Fig. [T] and 
m,f, in GeV. Both plots 



stau sector, M^-M L ^, see Eqs. ©, |[T0")|. The other MSSM 
parameters are given in Tabled] 



IV. SUMMARY AND CONCLUSIONS 



The CP-sensitive parts appear only in the spin-spin cor- 
relations, which can be probed by the subsequent neu- 
tralino decay 



ti + ii 



h -> x? 



(,, 



(50) 



for I = e, /x. The T-odd tau polarization normal to the 
plane spanned by the r and i\ momenta, can then be 
used to define a CP-odd tau polarization asymmetry. 
It is based on a triple product, which probes the CP 
phases of the trilinear scalar coupling parameter A T , the 
higgsino mass parameter /i, and the U(l) gaugino mass 
parameter Mi. 

We have analyzed the analytical and numerical depen- 
dence of the asymmetry on these parameters in detail. 
In particular, for nearly degenerate staus where the stau 
mixing is strong, the asymmetry obtains its maxima and 
can be larger than 70%. The normal tau polarization 
can thus be considered as an ideal CP observable to 
probe the CP phases in the stau and neutralino sector 
of the MSSM. 

Since the CP-sensitive parts appear only in the 
subsequent stau decay products the stau production 
process can be separated. Thus both, ILC, and LHC 
collider studies are possible. Concerning the kinematical 
dependence, the asymmetry is not Lorentz invariant, 
since it is based on a triple product. At the LHC, 
staus are produced with a distinct boost distribution. 
Evaluated in the laboratory frame, the resulting tau 
polarization asymmetries get typically reduced by a 
factor of two thirds, compared to the stau rest frame. 

We want to stress that a thorough experimental anal- 
ysis, addressing background processes, detector proper- 
ties, and event rate reconstruction efficiencies, will be 
needed in order to explore the measurability of CP phases 
in the stau sector at the LHC or ILC. We hope that our 
work motivates such a study. 
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We have analyzed the normal tau polarization and . ,. . ,„ , 

t ^-m 1 Appendix A: Momenta and spin vectors 

the corresponding CP asymmetry m the two- body decay 

chain of a stau 

For the stau decay f m — > T~Xi, we choose the coordi- 
fi r + X2- (49) nate frame in the laboratory (lab) system, such that the 
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momentum of decaying f points in the z-dircction. 

p>i = {E f ,0,0, |pH), (Al) 

p>{ = E T (l,sin0 T) O,cos0 T ), (A2) 
with the decay angle 8 T = <£(pf , p T ), and 

(m? -m| 9 ) 



« | Pr | W 2( ^_ _ | pf | c ; s0r) : 



(A3) 



in the limit m r — > 0. The momenta of the leptons from 
the subsequent neutralino decay \i — > I — > Xi^2 ©, 
can be parameterized by 

pf£ i = Ef 1 (1, sin 6*i cos 0i, sin 6\ sin 0i, cos #i ), (A4) 

p;jf = Ei 2 (l 7 sm02 cos</>2, sin# 2 sin0 2 , cos^), (A5) 
with the energies 

Et, 



2(EfO - |p£o|cos0 Dl )' 



(A6) 



1 1 

ra~ — 77i- 

Et, = ^ , , „ (A7) 



2(J5g- |p £ ~|cos6> D2 )' 

and the decay angles 0d 1 = <(Px?)P^i), 
<(p^,P£ 2 ), that is, 



cos Ud 1 



COS #D, 



(Pf - Pt) • gli 

Ipt-Pt| 

(gf - Pr - Pi i ) • Pl 2 

|p f - Pr - P£ X | 



the tau spin vectors by 
0, 



| S 2 X S|| 



= o, 



1 



Pr , | 7Pn 

Pr 



The spin vectors s°, a — 1, 2, 3, for the tau, and s~o, b = 
1,2,3, for the neutralino Xi> fulfil completeness relations 



Z — J 

a 

M v 



in 



(All) 
(A12) 



and they form orthonormal sets 

s a T -s c T = -S ac , s?.-p T = 0, (A13) 

4° -4? = 4° = °. ( A14 ) 

with p^ = p^/m. Note that the asymmetry A^ F , 
Eq. (|40|) . does not depend on the explicit form of the 
neutralino spin vectors, since they are summed in the 
amplitude squared, see Eq. ([3~T|) . using the completeness 
relation. 



Appendix B: Phase space 

The Lorentz invariant phase-space element for the stau 
decay chain, see Eqs. ^ - ((2]), can be decomposed into 
two-body phase-space elements [37| 

±Sfips(sf, p (l ,pe 2 ,Px°) = ?2^)2 d - SS - ps ( s ' : ' Pt > p x°J 

x ds^o dJz?ips(s^° ; p h , p~) ds^ dJ?ips(s£ p h , p^o). (Bl) 
The different contributions are 

*°) = T J Pt|2 2 sing T dg T , (B2) 



d£fips(s f ; p T ,Pj>o) 



dJZips{ Sn ;Pi x ,Pi) = 2( ^ )a m |T-m| dni ' (B3) 

^?P«(^^,Px?) = 2 (2^mf-i|o d " 2 ' (B4) 
with Sj = p 2 and d51 j = sin #j d0j d<^ . 

Appendix C: Density matrix formalism 

The coefficients of the stau decay matrix, Eq. (f2"5f , are 



(A8) 


D 


= y(Kil 2 + lCil 2 )(Pxr^) 


(A9) 




-S 2 JKe{<^ (&^ i )*}m^om T , 


We define 




= (+) y (l a ™| 2 - l 6 mil ) m r(Px? ' S r), 


X Pr \ 


yb 


= (+) y (Kml 2 ~ \Km\ 2 ) m x°(P-r ' 4?) 


X Pr|/ ' 


y-'flb 


= Y(i<-i 2 + i 6 ™i 2 )^-4?) m ^ 


(A10) 




+. 9 2 fHe{aL(&L)*} x 



(CI) 
(C2) 



(s a r - Pii o)(s b - r P T )-(s a T -4o)^o-Pr) 



- 5 2 Jm{a^(6^)*}[s^, p T , 4 ? , p^o]. (C4) 

The formulas are given for the decay of a negatively 
charged stau, f m — > t _ x°. The signs in parentheses hold 
for the charge conjugated decay t£, — ¥ r + Xi ■ 

Note that the terms proportional to m T in Eqs. (jCip . 
(|C2[) . and (|C4[) . are negligible at high particle energies 
E ^> m T , in particular can be neglected. 

The coefficients of the x? decay matrix, Eq. (|29l) . 
are H 



Di 



J/^| 2 Ko-to|), (C5) 
(Vl/ttl^^-ftO, (C6) 
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and the selectron decay factor is 

D 2 = ff 2 |/£| 2 K 2 -™x?)- ( C7 ) 

The signs in parentheses hold for the charge conjugated 
processes, that is Xi ~^ ^i^r m ^1- (EH- 

For the decay into a left slepton x_i — > ^t^L' Eq s - (ESJ), 
JC6|, and JC3) read jH 

Di = yl/^| 2 Ko-mf), (C8) 

V b Dl = { + ) 9 2 \fu\ 2 m^(s'i r p tl ), (C9) 

D 2 = g 2 \fh\ 2 H-m 2 -o), (CIO) 

respectively. The expressions for Eqs. (|31[) and (|32|) have 
to be changed accordingly. The sign in parenthesis in 
Eq. (|C9|) holds for the charge conjugated process xi ~> 



and the stau diagonalization matrix 1Z T , Eq. (IT2l . the 
Yukawa coupling 1^-, Eq. (p?!?|) . and the matrix [7, that 
diagonalizes the chargino matrix 0] , 

17* -M*± -V ] =diag(m-±,m-±). (D4) 

The stau decay width for the entire decay chain, Eqs. ([1]) 
- ([2]), is then given by 

r(r t44x?) = 




|A4| 2 dJz?ips(s f ; p T ,p tl ,p e2 ,p^o) (D5) 



= r(f) x BR(f -> tx?) x BR( X ° -> 

xBR(£^£ 2 Xi), (D6) 

(D7) 

with the phase-space element dJzfips, as given in the Ap- 
pendix El the amplitude squared 



Appendix D: Stau decay widths 



i\A(x»)\'\A(£)\'DD 1 D 2 , 



(D8) 



The partial decay width for the decay f m — > tx® in the 
stau rest frame is [31 1 



r(f m -> tx?) = -^r— !p A 



(Dl) 



with the decay function D given in Eqs. (|C1|) . and the 
approximation m r = 0. For the decay f m — > JA-xf the 
width is I3ll 



r(f m -> ^xf) = 



(to? — TO ~±) 2 
167TTO?' a 



" 2 |C,| 2 , (D2) 



with the stau-chargino-neutrino coupling 0, l3l| 

= -CCi)* Uji + Yt (K*T U f 2, (D3) 



obtained from Eqs. (f3TJ)) by summing the tau helicities A r , 
A' r . The neutralino branching ratios are given, for exam- 
ple, in Ref. [H], an d we assume BR(i -> I2X1) = 1- We 
use the narrow width approximation for the propagators 



|A(j)| 2 d s ,= 



mjTj 



(D9) 



which is justified for Tj/rrij <C 1, which holds in our case 
with Tj < 0(1 GcV). Note, however, that in principle 
the naive C(r/TO)-expectation of the error can easily re- 
ceive large off-shell corrections of an order of magnitude, 
and more, in particular at threshold, or due to interfer- 
ences with other resonant, or non-resonant processes j38| . 
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